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Feedback control schemes are a promising way to manipulate transport properties of driven col¬ 
loidal suspensions. In the present article we suggest a feedback scheme to enhance the collective 
transport of colloidal particles with repulsive interactions through a one-dimensional tilted wash¬ 
board potential. The control is modelled by a harmonic confining potential, mimicking an optical 
“trap”, with the center of this trap moving with the (instantaneous) mean particle position. Our 
theoretical analysis is based on the Smoluchowski equation combined with Dynamical Density Func¬ 
tional Theory (DDFT) for systems with hard-core or ultra-soft (Gaussian) interactions. For either 
type of interaction we find that the feedback control can lead to an enhancement of the mobility by 
several orders of magnitude relative to the uncontrolled case. The largest effects occur for interme¬ 
diate stiffness of the trap and large particle numbers. Moreover, in some regions of the parameter 
space the feedback control induces oscillations of the mean velocity. Finally, we show that the 
enhancement of mobility is robust against a small time delay in implementing the feedback control. 

PACS numbers: 05.60.Cd, 05.10.Gg, 02.30.Yy, 47.57.J- 


I. INTRODUCTION 

The manipulation of colloidal transport properties 
with feedback mechanisms has become an active topic 
of research in recent years. Examples include the im¬ 
provement of the net current in one-dimensional ratchet 
systems m, the transport of interacting particles in a 
tilted washboard under Pyragas control mis], the sorting 
of colloids in a micro-fluidic channel [7], and the adjust¬ 
ment of viscosity of a sheared colloidal suspension [5]. 
Further, feedback control has become an important con¬ 
cept in (bio-)particle trapping [SHH]. reaction-diffusion 
systems m, quantum transport [TBHT5] . laser dynamics 
m, and brain dynamics [MlHI]- An essential factor sup¬ 
porting the development of and theoretical research on 
feedback control strategies is the recent major progress 
of corresponding experimental techniques, among them 
the monitoring and steering of colloids [innHna and 
biomolecules [23], and preparation and non-destructive 
(weak) measurement of quantum states |M] . 

Within colloidal transport, most of the feedback stud¬ 
ies so far involve single particles mmm or dilute 
suspensions mm, ie., systems of non-interacting par¬ 
ticles. We note that, even in this idealized situation, 
feedback can induce effective interactions if the protocol 
involves system-averaged quantities m- For many real 
colloidal systems, however, direct interactions between 
the colloids stemming e.g., from excluded volume effects, 
charges on the particles’ surfaces, or (solvent-induced) 
depletion effects cannot be neglected. First studies of 
feedback control in presence of colloidal interactions in¬ 
dicate indeed complex dynamical scenarios. An example 
was considered in Refs. mm, where a Pyragas-type con¬ 
trol of colloidal transport in one dimension resulted in 
current reversal and oscillatory states. 

In the present work we explore the transport of in¬ 
teracting (repulsive) colloids in presence of a feedback- 


controlled harmonic “trap”. Indeed, trap-like devices ap¬ 
pear as a standard tool to implement feedback, both in 
experiments (see, e.g., [iiiMiiiKiaiisi) and in theory 
[25] . A prominent example is an optical laser tweezer act¬ 
ing on polarizable colloids. The corresponding trap po¬ 
tential can be modelled as a quadratic function in space 
[Ill 1291130]. 

In conventional applications the position of the center 
of the trap acting on the colloidal particle(s) is either 
constant in space, or it moves in an externally prescribed 
manner misiiiM]. In contrast to these situations (which 
are termed “open-loop” in control theory), we here con¬ 
sider a harmonic trap whose center coincides with the 
mean position of the particles. Thus, the trap potential 
depends on the particle’s position, yielding a feedback 
scheme. 

As a model system to demonstrate the principle of this 
feedback control we consider the paradigm example of 
colloids driven through a one-dimensional, spatially oscil¬ 
lating, tilted “washboard” potential with energy barriers 
much larger than the thermal energy. Already without 
feedback or any trap potential these (overdamped) sys¬ 
tems show interesting effects such as absolute negative 
mobility |33j or enhancement of diffusion at a certain 
“critical” force [HU [M]. Many single-particle transport 
properties in washboard potentials can be derived analyt¬ 
ically [55fl39j . Further, recent numerical studies indicate 
interesting interaction-induced transport phenomena, ex¬ 
amples being coherent motion of attractively interact¬ 
ing particles [3311311110], density excitations in Frenkel- 
Kontorova models mi, or single file diffusion [3S] . Given 
this background, one may expect that the interplay of 
external potentials, particle interactions, and feedback 
yields exciting additional effects. Our study shows that 
this is indeed the case. 

Specifically, we consider particles with either (in¬ 
finitely) hard or soft (Gaussian) repulsion, the latter de- 
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scribing polymeric particles in a coarse-grained fashion 
[42l |43]. The feedback control is implemented on the 
level of the Smoluchowski (overdamped Fokker-Planck) 
equation, in which the particle interactions are treated 
via Dynamical density functional theory (DDFT) [il] . 
Our numerical results demonstrate that the feedback- 
controlled trap in conjunction with particle interactions 
can lead to a drastic increase of the mobility by orders of 
magnitude. Loosely speaking, the particles “help to push 
each other over the energy barrier”. This phenomenon 
is accompanied by a freezing of the width of the den¬ 
sity distribution (thus opposing normal diffusion), and to 
time-periodic oscillations of the mean velocity not seen 
in the uncontrolled case. 

In the major part of our study we assume instanta¬ 
neous feedback. This is clearly an idealization given the 
fact that, in an experiment with feedback control, there 
is always one (or several) time delay(s) due to measure¬ 
ment, information processing, and implementation of the 
forcing |45) . However, the time delay of modern exper¬ 
imental feedback techniques for colloids [n [ini na is] 
is much smaller than the time scale of particle motion, 
justifying the approximation of instantaneous feedback. 
Still, to estimate the effects we also consider briefly the 
impact of time delay. 

The remainder of this paper is organized as follows. 
After the introduction of the theoretical background in 
Sec. |nj we discuss in Sec. [hT] the effect of feedback on a 
single particle. In the limit of vanishing washboard po¬ 
tential the transport can here be calculated analytically. 
The full problem is discussed in Sec. |IV[ where we present 
the main results. A conclusion is given in Sec. |V] 


II. MODEL 

We consider the motion of N interacting Brownian par¬ 
ticles in one dimension under the influence of the exter¬ 
nally imposed, tilted washboard potential 


Vext{x) = u{x) - X F , (1) 

where u{x) = uq sin^^nx/a). In Eq. Q, F denotes a con¬ 
stant driving force, and a and Uq denote wavelength and 
amplitude of the periodic potential u{x), respectively. On 
the particle level, the motion is described by the N cou¬ 
pled, overdamped Langevin equations 

-fXiit) = - u'{xi) -b F -b f™\xi, xn) 

+ ff^{xu.-.,XN) + ^2kBT-iUt) ( 2 ) 


—d/dxiJ 2 j^iV{xi,Xj) and is the force due to feed¬ 
back control, where DF stands for dynamic freezing. 
These forces are specified below. 

In the present study we rather describe the motion 
in terms of the space- and time-dependent one-particle 
density [44llT7] 


N 


Q{x,t) = x^{f)) 




( 3 ) 


where (...) denotes an average over all realisations of the 
random force The density is normalized according 

to f da; 0 (x, t) = N. The time evolution of g{x, t) is gov¬ 
erned by the extended Smoluchowski equation 


dtg{x,t) =^^dooxg{x,t) + -dx {g{x,t)dx[Ve^t{x) 

7 7 

-bVDF(a:, g) + Vi^t{x, £»)]) , (4) 


where the impact of particle interactions and of feedback 
control enters via the potentials Ijnt ^^nd EdF: respec¬ 
tively. Specifically, to treat the particle interactions we 
employ the concepts of Dynamical Density Functional 
Theory (DDFT) [44l [48l |49] . In this framework. 


Lint ( 3 ;, p) 


[^?] 

5g{x,t) ’ 


( 5 ) 


where Fint [£•] is the interaction part of an equilibrium free 
energy functional and 6/5g denotes a functional deriva¬ 
tive. Equation implicitly contains an adiabatic ap¬ 
proximation, i.e., the assumption that non-equilibrium 
correlations, at each time t, can be replaced by those of 
an equilibrium system with density g{x,t). 

We consider two types of interacting systems, that is, 
ultra-soft particles described by the Gaussian core model 
(GCM) and hard particles. The pair interaction potential 
according to the Gaussian core model reads 


VGCM(.Xi,Xj) = eexp 



( 6 ) 


This potential has been introduced as a coarse-grained 
(center-of-mass) potential between two fluctuating poly¬ 
mer chains, with the particle diameter a being propor¬ 
tional to the polymers’ radius of gyration [42l |43l [50]. 
To incorporate the GGM interactions into the dynami¬ 
cal equation Q, we employ the mean-field free energy 
functional 

\ Jdxdx' g{x,t) vggm{x- x') g{x',t). (7) 


for the position Xi(t) of the ith particle [36|. In Eq. ([^, 
the friction constant is denoted by 7 , Boltzmann’s con¬ 
stant by kB, the temperature by T, and the fi(t) are 
independent random numbers chosen from a Gaussian 
distribution with zero mean and unit variance. Eur- 
ther, /“* represents the force due to interaction be¬ 
tween particle i and other particles j, that is /“* = 


This functional has been proven to give a reliable descrip¬ 
tion of the equilibrium structure of the GGM, particu¬ 
larly at intermediate and high densities m- Gombining 
Eq. Q and ([^, we obtain 

£•) = y t) vgcm{x - x'). ( 8 ) 
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Hard particles with diameter cr are described by the 
interaction potential 


time. Solving Eq. (13) analytically with the initial con¬ 
dition p(x,t = 0 ) = S{x — xq) yields 


^hard ) — 


0 , for \xi — Xj\ > a 
oo , for \xi — Xj\ < a 


( 9 ) 


For one-dimensional systems of hard spheres there exists 
an exact free energy functional |52j 


T^hard 

inf 


If / fX+tTj2 '' 

[£>] = “2 y dx In ( 1 - y ^ dx'g{x',t) 


>^[Q{x+'^,t) + g{x-^,t)], ( 10 ) 


The free energy (10 1 corresponds to the one-dimensional 


limit of fundamental measure theory |53j . 

We now turn to the modelling of feedback control. To 
this end we use the potential 


VDF{x,g) =T]{x- {x)f. 


( 11 ) 


where 


{x){t) = ^ J dxxg{x,t) ( 12 ) 

is the time-dependent mean particle position. Thus, 
Eq. (Ill describes a moving harmonic trap centered 


around the mean position, resembling the potential 
seen by particles in moving optical traps [141 [29]. The 
strength of the harmonic confinement, g, is set to con¬ 
stant. Since Vdf depends on (x)(t) and, thus, on the 
dynamical state of the system, it corresponds to a true 
feedback control. This is different from an “open-loop 


controlled” moving trap, where (x) in Eq. (Ill would be 
replaced by a position moving with fixed velocity vq. 

We also note that the fact that our feedback control 
is coupled to a ensemble averaged quantity is in con¬ 
trast to other feedback mechanisms which are based on 
individual particle positions [15]. At the level of the 
Langevin equation ([^ our feedback control force reads 
/°^(xi, (x)) = 2g{xi - (x)(t)). Thus only depends 
on a single coordinate, x^, and on (x). This differs from 
other feedback control approaches where the feedback 
force itself depends on the number of particles 131231151. 


III. SINGLE-PARTICLE TRANSPORT 

To understand the basic properties of the effect of 
the feedback control Eq. (11) we first discuss the single¬ 


particle case {N = 1, J^mt = 0) without the periodic po¬ 
tential (ito = 0). In this case Eq. (|^ reduces to the one¬ 
dimensional Smoluchowski equation 


gdtg = k^Td^^g + - F)). 


(13) 


A main quantity characterizing the transport is the mean 
particle position (x), defined in Eq. (12 1 , as function of 


F 

(x)(t) = —t + Xo. 
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(14) 


Equation (141 shows that the mean particle position does 


not depend on the confinement strength g. This can also 
be seen by applying the coordinate transformation x' = 
X — vt to the Smoluchowski equation Eq. (13), setting 
V = F/g. With this transformation the term dx{gF) 
vanishes. Further, the force related to Vdf is invariant 
with respect to this transformation. Hence, the influences 


of F and g decouple. From Eq. (141 we calculate the 
mobility 


/i := lim 

t—^OO 
_ 1 

? 

7 


dt{x) 

F 


(15) 

(16) 


which only depends on the friction constant 7 . We will 
refer to this value of /i as the mobility of free motion. 

A further quantity of interest is the mean squared dis¬ 
placement 


w[t) = {{x-{x)f 


(17) 


To calculate w{t) we use (Vdf) = vx [cf. Eqs. (11) and 
(13)] which yields 


w{t) = ^ (1 - . 


2 g 


(18) 


For short times w{t) growths linearly with time, corre¬ 
sponding to diffusive behaviour. For long times diffusion 
is suppressed: w{t) approaches a limiting value deter¬ 
mined by g. Interestingly, a similar behaviour of wit) oc¬ 
curs in a model of feedback control of quantum transport 
m There, the fluctuations of the number of electrons 
tunneling through a quantum junction are suppressed 
with a feedback control force, which is linear in the fluc¬ 
tuation of the number of electrons. This corresponds 
to our harmonic confinement of the density fluctuation, 
and indeed, the two physically different situations are de- 
scribable by a formally identical Smoluchowski equation 

[ 52 ]. 

We now turn to the system in presence of the poten¬ 
tial u(x), defined below Eq. ([^. In the single particle 
case (iV = l, J^int = 0) Eq. ([^ then reduces to the Smolu¬ 
chowski equation 

gdtg = k^Tdxxg + 9xig(yuF + Kxt)) • (19) 


Without control {g = 0) Eq. (19) describes the thoroughly 
studied case of a Brownian particle in a washboard po¬ 
tential, where the mobility, as well as the long-time dif¬ 
fusion constant, are accessible analytically [36] 113 1^ . 
From that it is known that the mobility is very small if 
uq 3> A:bT and if the driving force F is smaller than the 
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FIG. 1. Single particle transport throngh a potential with 
wavelength a = 2.5a. (a-c) Density plots (black, left axis) and 
potential Vdf + u (grey, right axis). The time t is (a) 10® tb 
and (b,c) lO^^rs where the Brownian time tb = a’^^KksT). 
Part (d) shows the mean particle position with respect to 
time for different r;. For strong confinements an oscillatory 
behaviour emerges. 


so-called critical force Fc = uott/q, related to the diffu¬ 
sion maximum |37j (for F > the potential minima 
vanish). Otherwise the mobility is large, in particular it 
approaches I /7 for uq/F 0. The goal of our study is 
to enhance the mobility in the regime of deep wells. 


A. Numerical Results 


To explore the single particle transport for finite 77 and 
Vext ^ 0 we solve Eq. (19) numerically, choosing uq = 
ISfceT and F = 0.2Fc. As initial condition we choose 
the equilibrium (Boltzmann) distribution corresponding 
to the case F = 0 


q{x, 0) = exp + u{x))/kBT) jZ , (20) 

where Z is a normalisation constant. 

Figures [^a-c) show plots of the one-particle density 
g{x,t) for three values of r]. As expected for a trap, the 
width of the density distribution becomes the smaller the 
larger rj. Figure [^d) shows additionally the mean parti¬ 
cle position with respect to time. Interestingly, we find 
that at large values of 77 , oscillatory solutions emerge. 
At the corresponding values of 77 the conhnement is so 
strong that the particle is confined to a single well of the 
periodic potential, cf. Fig. [^c). 

We explain the occurrence of oscillations as follows. 
We take a view on the beginning of one step of an oscil¬ 
lation at time t = IO^tb for 77 = I/cbTct”^ [cf. Fig. [i|d )]. 
The potential Fdf + u at this time, shown in Fig. [Re) 
as grey shade, shows that the particle is localized at a 



FIG. 2. Velocity v{t) (a) and width w{t) (b) as function of 
time for different confinement strengths 77 for a = 2.5a. Os¬ 
cillatory solutions are characterized by sharp peaks, (c) Am¬ 
plitude Av and period T of the velocity as function of rj for 
a = 2.5a. Data points are shown only for values of rj where 
we find oscillations. The dashed line indicates the time I/tk 
where tk is Kramers’ rate. 


minimum of Fdf + u. As time progresses, the constant 
driving force causes the diffusion of the particle to the 
next minimum. This leads to a slow increase of the mean 
particle position (x). Then, the feedback control which 
moves with (x), lowers the energy barrier, and steadily 
accelerates the diffusion through the barrier. This leads 
to a fast motion until the particle arrives in the next well. 
The next oscillation then starts again with slow diffusion 
over the next barrier. The repeated cycle of motion con¬ 
sisting of slow and fast portions is particularly visible in 
the velocity 

v{t) = j^{x){t), ( 21 ) 

and the width w{t) which are plotted in Figs. [^a) and 
(b), respectively. Notice that 'w{t) oscillates around a 
constant value, reflecting that the width of the distribu¬ 
tion stays finite even at large times (“dynamic freezing”). 

We analyse the occurrence of these oscillations nu¬ 
merically in terms of period T and amplitude Av = 
(■Cmax — ■Cmin)/2 of velocity, shown in Fig. [^c). The val¬ 
ues Umax and Umin are the global maximum and minimum 
of v{t)\t>t,, respectively, where ti is a time after the dis¬ 
appearance of transients. From Fig. I^c) we find that 
oscillatory solutions occur in a range of intermediate 77 . 
In that range the amplitude Av increases with rj from 
nearly zero to large values. Furthermore, the period T of 
oscillations roughly coincides with the inverse Kramers 
rate, which is the relevant time scale for the slow barrier¬ 
crossing mentioned before. As we see in Fig. [Tl[d), the 
regime of pronounced oscillations partly coincides with 
a “speed up” of the motion. We quantify this “speed 
up” via an average mobility based on the time-averaged 
velocity 

1 

u=—y dtv{t) (22) 
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FIG. 3. (Color online) Single particle transport: Mobility 
in dependence on r] for a = 2.5a. The mobility is scaled 
by the mobility /ro of uncontrolled diffusion in a washboard 
potential. The feedback control can enhance the mobility by 
up to « 20%. Inset: Potential VuFix^t = 0) + u(x) for the 
three values of rj which are indicated by crosses in n{r]). 



FIG. 4. Single-particle transport in presence of the open- 
loop potential Vljpenioop [see Eq. (24|] with ?7 = lkB,Ta~^ and 
Vo = F/y at t = IOOtb and a = 2.5cr. Black line: one-particle 
density g{x,t). Grey shade: potential Popenioop + u. 


such that 

(23) 

Figure shows /r//io depending on r], where /tq « 1.2 x 
10“'^/7 is the mobility of the uncontrolled system {r] = 0) 
with the same external potential [36l[55]. For small rj, 
we find /r « The remaining deviation is presumably 
a numerical issue because, by definition, lim,,_>o/r = /tq. 
At intermediate values of t] the mobility shows a global 
maximum which lies above /tq. From comparison with 
Fig-i: c) we see that the maximum of lies in the 
range of rj where the oscillation periods of v(t) are about 
(in fact, somewhat smaller) than the inverse Kramers 
rate [361139j . Quantitatively, the maximal enhancement 
of mobility of « 20.4% is reached at 77 « 

For even larger values of 77 a sharp decrease of the mo¬ 
bility to zero is observed - the motion comes to a halt. 
To investigate this phenomenon we first note that the 
motion is always oscillatory (for these large 77 ) as long as 
there is transport at all [compare Figs. an d[^c)]. From 
the explanation of the oscillations above, we recall that 
the oscillation period is determined by the slow diffusion 
process over the energy barrier. The inset of Fig.j^shows 
the potential VbF( 2 ;, t = 0)+u{x) for three values of 77 . To 
ignite transport the particle must diffuse from the cen¬ 
tral valley at a: = 0 to the next valley at a; « 2.5cr. The 
larger 77 the larger the energy barrier. Thus, the larger 
77 the smaller the probability that the particle diffuses to 
the next valley, the longer the period of the oscillations, 
and the lower the mobility. For 77 = ik^T/a^ there is 
no motion at all in the time range of our calculations 
(t < IO^tb). 

Finally, we note that for single-particle transport the 
actual value of a is essentially arbitrary because a only 
determines the scales of time, density, and confinement 


strength, not the qualitative behaviour. 


B. Comparison with open-loop control 

To estimate the benefit of the feedback control scheme 
over the more established open-loop control, we briefly 
discuss the motion of a single particle under the potential 

Fopenloop — Xj {x Upt) , (24) 

where vq is a constant velocity of the trap. Choosing vq 
equal to the mean velocity v of the feedback controlled 
system, one observes the same general behaviour, but 
slight variations of oscillation frequency and amplitude. 
Large values vq > v lead, by construction, to a fast trans¬ 
port, but the particle is no longer located in the center 
of the trap. We can see this from Fig. [4] which shows 
the one-particle density for the velocity Vq = F/y [corre¬ 
sponding to free motion, see Eq. & and the effective 
trap generated by the potential I4penioop+u. In a real op¬ 
tical trap a large distance of the particle position to the 
center of the trap implies a large probability to escape 
pm. Hence, driving the particle too fast implies the 
risk of losing the particle completely. On the other hand, 
being too cautious and driving the particle too slowly 
is inefficient. Thus, the optimal velocity is difficult to 
predict in open-loop control. The feedback control auto¬ 
matically finds the optimal driving speed without taking 
the risk of losing the particle. Furthermore, the feedback 
control does not influence the direction of motion, it only 
enhances the absolute value of the mobility. 
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IV. MANY-PARTICLE TRANSPORT 

We now turn to interacting systems, as described by 
the SE 0 with Eq. 0 for ultra-soft particles and 
Eq. ( [To| ) for hard particles. There are now two relevant 
length scales, the wavelength of the periodic potential, 
a, and the particle diameter cr. Hence, the wavelength a 
is not just a scaling factor, as it was the case in single¬ 
particle transport. In addition, the number of particles 
N will play a role because the equations are now non¬ 
linear in g. In our numerical calculations, we set the 
ultra-soft particles’ interaction strength e appearing in 
Eq. (|^ to e = d/ce’E. The hard-particle interaction has 
no parameter besides a. The initial condition is set to 
the equilibrium density resulting at F = 0. 

In the following we study motion of clusters of inter¬ 
acting particles for various trap strengths 77 , numbers of 
particles N, and dimensionless wavelengths a/a. 


A. General behaviour 

The overall goal is to explore whether particle inter¬ 
actions enhance the efficiency of our feedback control in 
terms of the mobility. Before we start with the analy¬ 
sis of the mobility we want to give an impression of the 
general behaviour of our interacting systems. 

We begin our study of the three-dimensional parame¬ 
ter space (iV, rj, a) with small N and small rj. In Fig.j^c) 
we present a plot of the density profile and the potential 
VuF + u a.t r] = 0.01kBTa~‘^,N = 4. In fact, the density 
profiles shown in Fig.j^c) and Fig. j^a) are very similar. 
Similarities to the single-particle case vanish gradually 
if V or 77 (or both) are increased (at constant a), yield¬ 
ing larger values of the density in the trap. To describe 
the effect of these changes in density, we consider the 
effective potential Vint that one particle experiences due 
to the interaction with the other particles. The value 
of Vint at a position x increases with the corresponding 
densities g{x'),x' « x. Particularly large values of both, 
g{x) and Vint (a;), occur at the minima of Vdf + w- As a 
consequence, the potential Vdf + m -I- Vint, which governs 
the motion (together with the constant driving force), is 
characterized by smaller energy barriers than Vdf + u. 
Loosely speaking, Vint fills the valleys of Vdf + u [see 
Fig. I^b)]. For high densities, the hard particles form a 
“chain” and the ultra-soft particles form a cluster which 
is characterised by mutual overlap. In this dense situa¬ 
tion the contribution of Vi„t to Vdf + m -I- Vint can be¬ 
come so large that u becomes negligible. Thus, there are 
no hindering energy barriers any more. For both inter¬ 
acting systems, ultra-soft and hard particles, we actually 
find this case. Fig. [^a) shows the hard particle case as 
example. The potentials plotted in Fig. [^a) show that u 
in fact is a minor contribution to Vdf + w -I- Vnt ■ We con¬ 
tinue the discussion of parameter variations with focus 
on the mobility in Sec. |IVB[ 

Similar to the single-particle case, we find oscillatory 



30 40 50 60 70 80 90 x/a 


FIG. 5. (Color online) One-particle density q{x, t) on the 
left axis (black lines) and the potentials Vdf, u, and Vnt on 
the right axis. The parameters show (a) the case where u 
is suppressed, (b) the decrease of barrier height by interac¬ 
tion, and (c) the low- 77 /low-A regime. In detail, (a) shows 
A = 10 hard particles forming a chain in a narrow trap 
(77 = a = 2.5 ct) atf = 28.9rB, (b) A = 4 hard par¬ 

ticles forming a loose chain in an moderately narrow trap 
(77 = 0.7fcBTcr~'^,a = 2.5cr) at t = 985TB, and (c) A = 4 ultra- 
soft particles in a wide trap (77 = 10“^fcBTcr“^,a = 2.5a) at 
f=10®TB. 


solutions in the range of intermediate to large rj. For 
a representative system (A = 4 hard particles). Fig. 
summarises different characteristics of the oscillations in 
terms of width w{t), velocity v{t), and plots of the density 
for four times during one oscillation period. The oscilla¬ 
tion period is of the order of tb which is much shorter 
than the oscillation periods of several IO^tb we observed 
in the single-particle case [see Fig. ib)]. From Fig. [^we 
see that these oscillations are intimately related to con¬ 
figurational changes while the particle chain moves over 
a distance of about one wavelength a. Studying v{t) for 
different 77 , see Fig. we find that a couple of different 
oscillation patterns emerge. Moreover, the oscillations’ 
frequency rises with the mean of the velocity itself. This 
can be explained with the observation that the particles 
move one wavelength a during one period. Note that the 
maximal amplitude of oscillation neither coincides with 
largest 77 nor largest mean velocity. 


B. Mobility 


We now turn to the mobility, as a measure of the effi¬ 
ciency of feedback control. We define the mobility 77 in 
the same way as in the single-particle case via Eq. (231. 
Figures]^ and 1^ show 77 in dependence of 77 for ultra-soft 
and hard particles, respectively. 
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FIG. 6 . (a) One-particle density q(x, t) for N = 4 hard par¬ 
ticles subject to the tilted washboard potential I4xt [see part 
(b)] at rj = 3kBTa~^,a = 2.5a at four times. For clarity 
the curves are shifted vertically. Each time is labelled with 
a symbol, reappearing in parts (c) and (d) which show the 
width w{t), given by Eq. ( [T7| , and the velocity v{t), Eq. ( pT| ), 
over time, respectively. 



FIG. 7. (Golor online) Velocity v{t), Eq. ( [^ , for V = 10 
hard particles at a = 2.5cr for a broad range of rj. The ve¬ 
locity 3.7698cr/rB = E/y corresponds to the velocity of free 
motion under the force F. This value is reached at large rj. 
Oscillations occur for 77 > 0.05kBTa~^. 



FIG. 8 . (Golor online) Mobility p for ultra-soft particles in 
dependence of rj. Given that enough particles contribute, the 
mobility can rise to I/ 7 , the mobility of free motion. The 
thick line indicates the mobility in the uncontrolled case. Re¬ 
sults from Fig. are included with the notation A'’ = l. 



FIG. 9. (Color online) Mobility p for hard particles in de¬ 
pendence of 77 . The mobility reaches I/ 7 , the mobility of 
free motion, if enough particles contribute. The thick line in¬ 
dicates the mobility in the uncontrolled case. Results from 
Fig. are included with the notation N = l. 


For a > a we observe an extreme growth of fi with 77 
and N over several orders of magnitude for both parti¬ 
cle species. We explain this behaviour with the corre¬ 
sponding decrease of the height of the energy barriers in 
Vdf + w + Mnt, which results in a larger diffusion rate 
and a faster transport. The same effect was observed in 
a study of the transport of super-paramagnetic colloids 
[5Ii| . For certain 77 and /i increases even up to the max¬ 
imal possible value 77 = I/ 7 , the mobility of free motion. 
An example for this large mobility is the case of = 10 
hard particles at 77 = IO/cbTct”^, shown in Fig. [^a). In 
this case, there are no hindering energy barriers (as we 
have analyzed in Sec. IV AI which then results in the high 
mobility. To achieve this suppression of u the well cre¬ 
ated by the trap potential must be very deep, i.e. SOOfeeT 


for the exemplary case shown in Fig.[^a). This value ex¬ 
ceeds those in typical experiments with light fields [HIST]. 
However, the transport of A^ = 40 ultra-soft particles at 
the mobility /i = I /7 at a = Scr needs a trap which is 
only AOk^T deep. 

Further, we see from Figs. that our feedback con¬ 
trol does not lead to a significant speed up for a = a. By 
analysing the potential landscape for a = <j for a series of 
77 and N (not shown) we find that the effective potential 
Vint develops peaks between the minima of Fdf + w. This 
means that the effective potential barrier encountered by 
a moving particle increases when 77 or A^ is enlarged. 
This is in contrast to the case a > a where the peaks of 
Vint are found at the minima of Vdf + u [see Fig. [^b)]. 
Our interpretation for the case a = a therefore is that 






























































































































the particles “pin” each other to the potential minima of 
u{x). 

We now consider the behaviour of for large 77 

(Figs. EE- For small N we observe a breakdown of 
motion, similar to the one observed in the single-particle 
case (see Fig. [^. However, this breakdown is shifted to¬ 
wards larger value s of rj . We recall that an increase of 
N at fixed 77 (Sec. IV A) leads to a decrease of the bar¬ 
riers of the potential Vdf + m -I- Vnt- This enhances the 
mobility (relative to that at V = 1 ) in the first place. 
However, upon increase of 77 (at fixed N) there can be 
a situation where the diffusion rate is not sufficient any 
more to populate the next local minimum of the potential 
Vdf + w. This is where transport breaks down. The com¬ 
bination of these two effects leads to the observed shift of 
the breakdown of mobility. Upon further increase of N 
and 77 , there comes a point where the large energy scales 
of Vdf and Vint suppress any influence from u. There¬ 
fore, we expect that the transport for high N exists for 
arbitrarily large rj. 

In Fig. we see that the increase of a at constant 77 
and N leads to an enhancement of mobility (as long as 
there is transport at all). This can be explained with 
the potential Vdf + u, whose valleys become broader the 
larger a. In a broader valley more particles accumulate 
which strengthens the role of interaction for the barrier 
crossing. However, this effect is limited by N: The parti¬ 
cle number must be large enough to fill at least one valley 
with particles, otherwise the transport breaks down. 


C. Time delay 


In a realistic set up with feedback control, a finite time 
is required to perform the measurement required to de¬ 
fine the control (In the present case, this measurement 
process concerns the average particle position). Hence, 
there is a certain time delay TDeiay To explore the sensi¬ 
tivity of our results towards TDeiay we change the control 
potential given in Eq. ( |11[ ) into the expression 

VDF’'^(a:, g)=r]{x- {x){t - TDeiay))^ ■ (25) 


We now consider two special cases involving hard parti¬ 
cles, where the non-delayed feedback control leads to a 
particularly high mobility (see Fig.[^. Numerical results 
are shown in Fig. The delay causes a pronounced 
decrease of mobility which appears to be linear in TDeiay 
for small delay times. Realistically, feedback mechanisms 
can be implemented at the time scale of 10 ms m^m 
whereas tb, the timescale of Brownian motion, is for 77 m- 
sized particles in the order of 300 ms [1] or larger [58l [59] . 
Hence, we expect that the ratio TDeiay /^b is rather small, 
that is, of the order 10“^. For such situations, our results 
in Fig. IT predict only a small decrease of 77 relative to 
the non-delayed case. However, even for large delays the 
mobility only decreases about one order of magnitude. 
This implies that even the time delayed feedback con- 



TOelay/TB 


FIG. 10. Effect of time delay on mobility 77 for hard particles. 
Even for very long delays the mobility is only reduced by one 
order of magnitude. 


trol can enhance the mobility by more than two orders 
of magnitude with respect to the uncontrolled case. 


V. CONCLUSION 

Inspired by the physics of a moving optical tweezer, 
we have proposed a feedback control strategy for the col¬ 
lective transport of interacting colloids through a corru¬ 
gated potential landscape. Our main goal was the theo¬ 
retical demonstration of the working principle for a well- 
defined model system. To this end we have considered 
the one-dimensional, overdamped motion of colloids with 
either hard or soft repulsive interactions in a tilted wash¬ 
board potential. The feedback control enters into the 
(Smoluchowski-like) equation of motion via a harmonic 
potential centered at the mean particle position. Thus, 
contrary to other studies [23 [23 ) the present feedback 
control cannot induce motion on its own. 

The main result of our study is that the interplay of the 
feedback control, on one hand, and particle interactions, 
on the other hand, can generate a drastic increase of the 
average mobility by several orders of magnitude relative 
to the uncontrolled, single-particle reference case. The 
largest mobilities occur for rather stiff traps and high 
densities (i.e., large N) inside the trap, yielding chain- or 
cluster like packages of colloids. Here, the mobility rises 
up to its limiting value defined by the mobility of a freely 
moving, overdamped particle. Interestingly, this giant 
increase does not occur for a single particle under the 
same feedback control. This shows that the observed mo¬ 
bility enhancement is indeed an interaction effect. The 
enhancement can be explained by the fact that, in pres¬ 
ence of particle interactions, these dominate the effective 
“field” acting on an individual particle, while the impact 
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of the external potential barriers vanishes. Thus, parti¬ 
cles “help each other” to overcome the external barriers. 
Another new feature is the emergence of oscillatory be¬ 
havior of the mean velocity (and the width of the density 
distribution) due to the feedback-controlled trap. The 
latter effect occurs for both, single and interacting par¬ 
ticles, with the period of oscillations being close to the 
inverse of Kramers’ escape rate. 

From an application point of view it is interesting that, 
due to its coupling to the mean position, the feedback- 
controlled trap implies a small risk to “lose” particles. 
Indeed, the width of the distribution stays constant on 
time-average, reflecting a “dynamical freezing”. This 
is different from externally moved, “open-loop” traps, 
where an inappropriate choice of the trap velocity eas¬ 
ily lead to a broadening of the density distribution, and 
thus, a sprea ding o f particles out of the trap (see dis¬ 
cussion in Sec III B). Another experimentally relevant is¬ 
sue concerns the impact of time delay(s). Here we have 
shown that time delay does indeed reduce the mobility, 
similar to what has been observed in ratchet systems [3] . 
However, for realistic time delays the remaining mobility 
is still enhanced by two orders of magnitude. 

Concerning the methodology, we note that the DDFT 
scheme employed here implies an “adiabatic” approxima¬ 
tion of the time-dependent two-particle correlations. It is 
now well established gHUSOlISI] that this approximation 
may generate artefacts especially for densely packed par¬ 


ticles, e.g., during the expansion of a cluster. Since we 
are mainly focusing on steady transport conditions we ex¬ 
pect our results to be at least qualitatively right. Still, it 
would be very interesting and important to test our pre¬ 
dictions against explicit Brownian Dynamic simulations 
of the corresponding (overdamped) Langevin equation, 
Eq. ([2|. 

Finally, we would like to point out that the con¬ 
cept behind dynamical freezing is not restricted to one¬ 
dimensional washboard potentials. Indeed, the present 
feedback control can easily be formulated in two or three 
spatial dimensions. Further, the external potential hin¬ 
dering the motion does not have to be static or even peri¬ 
odic which enriches possible applications. An interesting 
question is how well the present control strategy works 
for other types of colloidal interactions, particularly at¬ 
tractive ones. Another open question concerns the impli¬ 
cations for the (non-equilibrium) thermodynamics of the 
system, an area which currently receives much attention 
[4l [62l[64| . Work in these directions is in progress. 
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